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^  Graphical  Techniques  for  Maintenance  Planning 

Items  subject  to  stochastic  failure  are  candidates 
for  preventive  maintenance  where  often  an  item  is  replaced 
or  restored  to  a  good-as-new  condition.  For  single  com¬ 
ponent  items,  the  optimal  interval  for  this  type  of  pre¬ 
ventive  maintenance  might  be  determined  via  some  mathe¬ 
matical  model  of  replacement.  Usually,  these  models 
require  knowledge  of  the  failure  distribution  as  well  as 
the  costs  of  preventive  and  corrective  maintenance.  For 
most  problems  there  is  considerable  uncertainty  as  to  the 
appropriateness  of  a  failure  model.  Bergman  and  others 
have  espoused  graphical  tp'ihniques  for  the  Age  Replacement 
model  which  do  not  require  specification  of  the  failure 
model.'  We'^ropos^  several . graphical  techniques  for  two 
other  replacement  models.  Block  Replacement  and  one  w&  calleJl 
Blind  Replacement,  cs" 

These  techniques  use  historical  data  (a  moderate 
cunount)  to  estimate  the  optimal  replacement  interval.  In 
the  absence  of  such  data  the  techniques  here  do  not  apply. 
With  a  view  toward  small-sample  size  statistical  inference, 
we  investigate  the  bias  in  graphical  estimates  of  the 
optimal  replacement  interval,  t*,  for  Age,  Block,  and  Blind 
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Replacement.  We  find  that,  for  small  samples,  techniques 
using  the  total  time  on  test  statistic  may  have  positive 
bias.  We  propose  a  survivor  time  on  test  statistic  which 
results  in  estimates  that  may  have  a  negative  bias.  We 
recommend  composite  estimators  in  some  cases  to  mitigate 
bias  when  sample  sizes  are  small. 

These  graphical  techniques  provide  only  point 
estimates  at  t*.  Accordingly,  we  propose  two  confidence 
bounds  techniques  using  a  bootstrap  resampling  approach. 
For  Age  and  Blind  Replacement,  we  show  that  these  tech¬ 
niques  in  many  cases  produce  useful  confidence  bounds 
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CHAPTER  I 


INTRODUCTION  AND  SUMMARY 

Introduction  to  Planned  Maintenance 
Maintenance  strategies  are  sets  of  actions  to  be 
taken  periodically  on  a  repairable  system  in  order  to  sus¬ 
tain  its  performance.  An  optimal  strategy  results  in 
system  performance  at  its  best,  as  measured  by  a  specified 
criterion.  System  performance  can  be  classified  into  dis¬ 
crete  categories  of  conditions,  such  as  a  two-state  case 
of  good  and  failed  or  a  multi-state  case  with  conditions 
between  these  extremes.  For  simple  systems,  that  is 
single-component  systems  which  either  are  in  good  condi¬ 
tion  or  have  failed,  we  can  consider  four  basic  mainte¬ 
nance  strategies.  Two  of  these  strategies  apply  in  the 
situation  when  the  system's  condition  is  known  at  every 
instant,  while  the  other  two  apply  when  the  condition  is 
not  known  except  by  inspection. 

The  first  basic  strategy,  used  when  the  system's 
condition  is  known,  is  that  of  replacing  (or  restoring  to 
a  "good-as-new"  condition)  the  item  at  failure  or  at 
age  t,  whichever  occurs  first.  This  policy  is  known  as 
Age  Replacement  and  requires  knowledge  of  both  the  item's 
age  and  its  condition.  The  second  strategy  uses 


information  on  the  item's  condition,  not  its  age;  this 
policy  is  to  replace  the  item  at  failure  and  at  fixed 
points  in  time.  Such  a  policy  is  known  as  Block  Replace¬ 
ment  ;  it  avoids  the  expense  of  recording  failure  times  and 
tracking  age.  However,  it  incurs  the  additional  cost  of 
possibly  replacing  relatively  new  items. 

In  the  above  two  strategies,  the  only  uncertainty 
is  that  of  an  item's  "lifetime,"  i.e.  the  age  when  failure 
occurs.  On  the  other  hand,  if  an  item's  condition  (either 
good  or  failed)  is  unknown  except  by  periodic  inspection, 
then  there  are  two  kinds  of  uncertainty;  that  of  the 
item's  lifetime  and  that  of  the  item's  condition  between 
inspections.  The  third  and  fourth  basic  maintenance 
strategies  address  both  of  these  uncertainties.  The  third 
strategy,  which  we  call  Maintenance  Inspection,  is  to 
inspect  the  item  periodically  every  t  units  of  time  and 
to  replace  only  if  there  is  a  discovered  failure.  The 
fourth  strategy,  which  we  refer  to  as  Blind  Replacement, 
is  to  replace  the  item  regardless  of  its  condition  every 
t  units  of  time.  This  fourth  strategy  is  prudent  if  the 
cost  of  replacement  is  less  than  the  cost  of  inspection 
(as  would  be  the  case  if  an  inspection  were  destructive  to 
the  item) •  Note  that  under  both  Maintenance  Inspection 
and  Blind  Replacement,  a  failed  unit  will  not  be  replaced 
at  the  time  of  failure.  These  strategies  are  to  be  used 
only  under  circumstances  in  which  the  failure  of  the  item 
is  not  fatal  to  the  system. 


In  this  chapter,  we  provide  a  synopsis  of  the 
literature  on  three  of  these  four  strategies  and  summarize 
our  findings  and  conclusions.  Since  Maintenance  Inspec¬ 
tion  models  are  complex  {cf.  Brender  (1963)}  and  do  not 
appear  amenable  to  graphical  solution,  we  do  not  address 
them. 

All  of  the  three  strategies  that  we  address  have  a 
common  structure  known  as  a  renewal  process.  In  a  renewal 
process,  the  instant  in  time,  at  which  the  item  is  made 
"good-as-new"  by  either  repair  or  replacement,  is  a  renewal 
point.  A  sequence  of  renewal  points  through  time  is  a 
renewal  process  as  shown  in  Figure  1. 


renewal  renewal  renewal  renewal 

time  tine  time  time 

Fig.  1.  A  General  Renewal  Process 


The  times  between  renewal  points,  for  i=l,  ...  which  are 

random  variables,  can  be  summed  to  obtain  the  time  of  the 
th  ” 

n  renewal  point,  S  =  EL..  Moreover,  there  is  usually 

"  i=l  ^ 

a  cost  (of  repair,  replacement,  inspection,  unavailability, 

etc.),  V^,  incurred  during  S^]  the  i^*'  renewal 

interval.  This  cost  is  also  a  random  variable.  The 

average  cost  per  unit  time  during  the  interval  (0,  S  ] 

n  n  ” 

equals  the  ratio  E  V./  E  L..  Assuming  that,  in  the  two 

i=l  ^  i=l  ^ 

sequences  of  random  variables  {V^,  ...,  V^,  ...}  and 
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{Lj^,  L^,  (V^,  L^)  are  independent  and  iden¬ 
tically  distributed,  Ross  (1983:78-79)  proves  that  the 
long  run,  average  cost  per  unit  time  is 


where  e{v}  =  e{Vj^}  =  ...  =  eCv^^)  =  ...,  and  e{l}  =  EtL^^) 

=  ...  =  e{l^}  =  ...  .  The  solutions  techniques  for  the 
probability  models  which  we  address  all  seek  the  value  of 
the  decision  variable  which  minimizes  this  long  run, 
average  cost  per  unit  time. 

Age  Replacement  Model 

Barlow  and  Proschan  (1965)  discuss  Age  Replacement 
using  long  run,  average  cost  per  unit  time  as  the  cri¬ 
terion.  They  presume  a  failure  results  in  a  cost,  Cj^, 
and  that  each  replacement  incurs  a  cost,  C^t  where  either 
a  failure  or  a  replacement  results  in  a  renewal  of  the 
system.  They  show  that  the  long  run,  average  cost,  C(t), 
for  a  policy  of  replacement  at  failure  or  at  age  t  can  be 
expressed  as 

C  F(t)  +  C-  R(t) 

C(t)  =  — - ^ »  (1.2) 

q/  R(x)dx 


where  F(t),  the  cumulative  distribution  function  (cdf ) , 
is  the  probability  that  a  failure  will  occur  prior  to  or 
at  age  t,  and  R(t)  is  the  probability  that  the  item  will 
survive  beyond  age  t  {i.e.  R(t)  *  l-F(t)}.  The  numerator 
of  equation  (1.2)  is  the  expected  cost  per  renewal  inter¬ 
val,  E{V},  while  the  denominator  represents  the  .expected 
length  of  a  renewal  interval,  E{L),  given  that  lifetimes 
are  truncated  at  age  t,  i.e. 

t  t 

E{min(L,t)}  =  /  xdF(x)  +  t  /  dF(x)  =  f  R(x)dx  . 

0  t  0 

The  quantity  f  R(x)dx  =  y(t)  has  a  crucial  role  in  Age 
0 

Replacement.  It  is  the  expected  operating  time  of  an 
item  which  is  withdrawn  from  service  upon  attaining  age  t 
Clearly,  y  (<»)  is  the  mean  life  (or  mean  time  between 
failures  given  that  there  is  no  preventive  replacement) 
while  y(t)  is  the  "conditional  mean  life”  (or  mean  time 
between  maintenance) . 

Solution  to  Age  Replacement 
When  F(t)  is  Specified 

Since  R(t)  =  l-F(t),  equation  (1.2)  can  be  alge¬ 


braically  rewritten  as 
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Although  the  criterion,  C(t),  is  recast  in  equation  (1.3) 
as  the  long  run,  average  standardized  cost  per  unit  time, 
a  graphical  clarity  can  be  obtained  as  shown  in  Figure  2. 

If  we  consider  the  curve  {x(t),  y(t))  =  {p(t),  F(t)}  for 
0  £  t  00,  the  slope,  c,  of  any  line  y  =  cx-C^/ (Cj^-C^) 
from  the  point  {0,  -C2/ (Cj^-C2)  )  to  {x(t),  y(t)}  equals 
C(t).  The  tangent  line,  as  drawn  in  Figure  2,  is  the  line 
with  minimum  slope  and  thus  represents  the  optimal  solu¬ 
tion.  The  parameter  value  at  the  point  of  tangency,  t*, 
is  the  Age  Replacement  interval  which  achieves  a  minimum 
for  C(t).  Therefore,  t*  is  an  optimal  interval  with 
C(t*)  equal  to  c* ,  the  slope  of  the  tangent  line.  Note 
that  there  may  be  more  than  one  point  of  tangency,  in  which 
case  there  are  multiple  solutions.  If  p  (00)  <  there 
will  always  be  at  least  one  point  of  tangency  as  can  be 
seen  from  the  geometry  in  Figure  2.  A  nice  feature  of  this 
graphical  approach  is  that  it  avoids  regularity  conditions 
for  unique  optimum  solutions  such  as  F(t)  belonging  to  the 
Increasing  Failure  Rate  (IFR)  class  of  distributions.  We 
remark  that  the  curve  {x(t),  y(t))  is  increasing  (that  is 
to  say  nondecreasing)  in  t  because  both  F(t)  and  ij(t)  are 
increasing  functions  in  t. 

Solution  to  Age  Replacement 
When  F(t)  is  Estimated 

In  order  to  use  equation  (1.3)  and  its  graphi  al 
representation  for  determining  t*,  the  failure  distribu¬ 
tion,  F(t),  must  be  specified.  Often  in  practici',  F(t)  is 


not  known,  and  neither  is  it  reasonable  to  assume  a  life¬ 
time  distribution.  Arunkumar  (1972)  addresses  this  prob¬ 
lem  using  ordered  failure  times  to  arrive  at  an  empirical 
failure  distribution  and  at  an  analytical  solution  for  t* 
Bergman  (1977a)  develops  a  graphical  technique  for  esti¬ 
mating  t*  using  the  "total  time  on  test  statistic." 

Bergman's  technique  supposes  that  there  are  n 
observed  lifetimes  {t^^,  •..»  t^}  of  an  item.  These 

are  "free-flowing"  (unrestricted)  life  lengths  from  a 
complete  life  test  obtained  from  the  laboratory  instead 
of  from  truncated  field  data.  These  n  lifetimes  can  be 
ordered  as  ^ ^  (i)  1  t  (2)  1  •  •  •  1  ^  (n)  ^ "total 
time  on  test,"  through  the  i^^  failure  time,  t^^ 

can  be  calculated  as 


where  tQ  =  0.  The  ratio  T(t  ) /T  (t^^j )  =  U^  is  the  scaled 
total  time  on  test  at  age  and  a  plot  of  versus 

i/n  is  known  as  a  total  time  on  test  (TTT)  plot.  As  we 
discuss  below,  Bergman's  idea  is  to  change  the  scale  of 
the  denominator  of  equation  (1.3)  so  as  to  take  advantage 
of  the  geometric  properties  of  the  total  time  on  test  plot. 
He  then  approximates  this  scaled  denominator  by  and 
estimates  the  failure  distribution,  in  the 

numerator  of  equation  (1.3)  by  i/n. 

Bergman's  technique  is  illustrated  in  Figure  3  for 
a  sample  of  size  four  wherein  we  plot  the  point  '^2^ 
on  the  horizontal  axis  and  construct  a  tangent  to  the  TTT 
plot  passing  through  this  point.  We  remark  that  Figure  3 
is  a  reflection  of  Figure  2  rescaled  about  the  line  y  =  x. 

We  present  Bergman's  technique  as  it  is  in  Figure  3  because 
others  have  used  this  same  graphic  portrayal  [cf.  Barlow 
(1978)].  For  reasons  discussed  below,  we  prefer  the 
graphical  representation  in  Figure  2 . 

Bergman's  technique  involves  several  approxima¬ 
tions.  First,  he  estimates  the  theoretical  failure  dis¬ 


tribution  by  an  empirical  distribution,  F^.  Since 

lim  F  =  F,  a  finite  number  of  observations  results  in  an 
n 

n-^-oo 

approximation  of  F  that  improves  as  the  number  of  obser¬ 
vations  is  increased.  With  ample  observations,  then. 


equation  (1.3)  can  be  approximated  as 


c:^  ^  ''n“' 


C(t)  '  C^(t)  = 


/  R  (x)dx 
n 


(1.4] 


Bergman  (1977a)  states  that  an  estimate  of  t*,  the  mini- 
mizer  of  equation  (1.4) ,  may  be  found  aunong  the  n  ordered 
lifetimes  {t,,.  <  ...  <  t.  "thus  to  estimate  the 

optimal  age  replacement  interval  it  is  enough  to  find  the 
index  j  for  which 


^n<Nj))  = 


Cj^-C2  ^n‘^(j)^ 


(1.5 


R^(x)dx 


is  a  minimum."  Finally,  Bergman  estimates  the  empirical 
failure  distribution  using  the  plotting  position  "i/n"  and 
the  conditional  mean  life  using  the  total  time  on  test 
statistic,  divided  by  n.  A  plotting  position  is 

a  representation  of  the  empirical  cumulative  probability, 

F  (t,.,),  of  the  i*^^  ordered  observation,  t,...  We  dis- 
cuss  the  plotting  position  of  i/n  in  greater  detail  later. 
Thus,  using  these  two  estimates,  equation  (1.5)  is  further 
approximated  as 


(1.6) 


(1.7) 


Since 


In  lieu  of  minimizing  equation  (1.7)  Bergman  maximizes  its 


reciprocal,  i.e.  he  maximizes  the  ratio 


since  the  term  n/T(tj^j)  is  constant.  The  geometric  idea 
is  that  this  ratio,  which  is  the  slope  of  any  line  passing 
through  the  point  {-C2/ (€^^-€2)  »0 )  and  the  TTT  plot,  is 
proportional  to  Cj^(t^j)  in  order  to  minimize 

Cn(t(i)),  we  need  only  maximize  this  slope.  The  maximum 
slope  is  attained  by  a  line  tangent  to  the  TTT  plot  pass¬ 
ing  through  the  point  { -C2/ (C ^^-€2 )  ,  0  }  as  in  Figure  3. 

In  Figure  3,  the  numerator  i  in  the  value  i/n  correspond¬ 
ing  to  the  point  of  tangency  is  the  index  of  ,  an  esti¬ 

mate  of  the  optimal  replacement  interval,  t*  (in  Figure  3 
this  estimate  is  the  failure  time  of  the  3rd  failure) . 

Bergman's  technique  is  a  significant  contribution 
to  the  area  of  maintenance  planning.  With  this  method, 
maintenance  analys'  s  can  avoid  the  necessity  of  goodness- 
of-fit  testing  and  the  inherent  risk  of  incorrectly  speci¬ 
fying  a  theoretical  failure  distribution.  Although  our 


!  t 


Fig.  3.  Bergman's  Graphical  Technique 


graphical  approach,  in  Figure  2,  might  appear  at  first  to 
differ  from  Bergman's  approach,  each  is  a  reflection  image 
of  the  other  (with  different  scaling)  which  leads  to  the 
same  result.  Thus,  our  rationale,  embodied  in  Figure  2, 
is  only  a  refinement  of  Bergman's  graphical  argument. 

From  a  practical  point  of  view,  Bergman's  method 
is  appealing  because  it  lends  itself  to  sensitivity  analy¬ 
sis  for  cost  values.  Since  management  may  be  more  com¬ 
fortable  with  a  range  of  cost  values  rather  than  point 
values,  tangent  lines  drawn  through  the  end  points  of  a 
range  of  costs  will  display  cost  sensitivity  in  terms  of  a 
range  of  the  estimated  optimal  Age  Replacement  interval. 


Plotting  Position  and  the  Total 
Ti—  on  To»t  Statistic 

Bergnan'a  technique  uses  the  plotting  position 

i/n,  knoim  as  the  California  position,  to  estinate  the 

failure  distribution  at  California  position  is 

perhaps  the  most  frequently  encountered  representation  of 

F_,  and  as  we  show  below,  it  is  intrinsic  to  the  total  tisie 
n 

on  test  statistic.  There  are  other  plotting  positions 
such  as  i/Cntl),  the  mean  position;  (i-l)/(n-l),  the  modal 
position;  and  (i-l/2)/n,  known  as  Hazen's  position. 

Weibull  ( 1939a, b)  reconanended  the  mean  position  since,  for 
the  i^^  order  statistic  (a  random  variable),  the 

expectation  of  the  cumulative  probability  at  the  i^^  order 
statistic,  E{F(X^j)},  is  i/(ntl).  This  result  holds 
regardless  of  the  underlying  distribution  generating 
ordered  observations.  Harter  (1984)  provided  a  comprehen¬ 
sive  rev'iew  of  the  plotting  position  literature  and  noted 
that 

. . .  much  of  the  disagreement  and  confusion  as  to  the 
choice  of  plotting  positions  is  due  to  the  fact  that 
the  cdf  at  the  expected  value  of  the  ith  order  sta¬ 
tistic  is  not  equal  to  the  expected  value  of  the  cdf 
at  the  ith  order  statistic,  i.e. 

riE{X(i)  }  ft  E{F(X(i))  } 

except  in  the  case  of  the  uniform  distribution 
F(x)  =  X,  0  <  X  <  1. 

In  Chapter  II,  we  investigate  the  question  of  bias 
in  Bergman's  technique  via  Monte  Carlo  simulation  and  con¬ 
clude  that,  for  some  common  distributions,  it  produces 
biased  estimates  of  t*.  There  may  be  several  reasons  for 
this  bias.  For  example,  Cunnanc  (1978)  shows  that  the 


California  plotting  position  eimong  others  is  biased  for 
several  common  probability  plots.  Such  a  bias  would 
affect  Bergman's  technique  both  in  its  estimate  of  F(t) 
and  of  y (t) . 

Any  bias  in  y (t)  would  stem  from  the  total  time 
of  test  statistic  which  is  directly  related  to  the  empiri¬ 
cal  failure  distribution,  F  (t) .  The  standard  convention 

n 

is  to  represent  F^(t)  as  a  right  continuous  step  function 
(as  we  show  in  Figure  4  for  the  uniform  distribution  and 
its  California  position  representation) .  The  plotting 
position  determines  the  step  size  at  t  =  tj^j  and  the 
starting  value  F^(t=0).  In  Figure  4,  we  depict  the  steps 
occurring  at  =  i/ (n+1)  for  ease  of  illustration, 

although  in  practice  this  would  be  an  unlikely  realiza¬ 
tion.  The  standard  empirical  survival  distribution, 

R  (t) ,  on  the  other  hand,  equals  1-F  (t)  as  shown  in 
n  <  n 

Figure  5.  Regardless  of  the  shape  of  the  distribution,  an 

estimate  of  y(t,..),  which  is  the  area  under  R  (t,..), 
vij  n  (i( 

can  be  calculated  for  the  California  plotting  position  as 

n 
n 

n 


(1) 


n-1 


*'(2)~*'(1) 


n-i+1 


^(i)"’'(i-l) 


where  tQ  =  0.  Note  that  the  numerator  is  the  total  time  on 
test  statistic,  T(t,.,).  It  follows,  then,  that  any  bias 


L5 


in  the  empirical  failure  distribution  due  to  the  Cali¬ 
fornia  plotting  position  will  also  bias  the  estimate  of 
mean  life,  )j(t),  through  the  total  time  on  test  statistic. 

As  an  alternative  to  the  total  time  on  test  sta¬ 
tistic,  we  propose  another  statistic  which  we  call  the 
"survivor  time  on  test  statistic,"  S(tj^j).  This  sta¬ 
tistic  is  based  on  the  time  on  test  accumulated  by  sur¬ 
viving  items  which  live  beyond  the  age  .  It  is  com¬ 

puted,  at  the  i^^  failure  time, 

where  tQ  =  0.  In  Chapter  II,  we  show  that  a  composite 
estimator  of  t*  using  both  and  S(t^^)  may  produce 

less  bias  under  some  circumstances  than  Bergman's  esti¬ 
mator. 

Block  Replacement  Model 

Barlow  and  Proschan  (1965)  also  discuss  a  proba¬ 
bility  model  for  the  second  maintenance  strategy  where  an 
item  is  replaced  at  failure  as  well  as  at  fixed  calendar 
time  intervals  of  length  t  regardless  of  age.  Their  objec¬ 
tive  is  to  minimize  long  run,  average  cost  per  unit  time 
where,  as  in  their  Age  Replacement  model,  Cj^  denotes  the 
cost  of  a  failure  and  represents  the  cost  to  replace 
the  item,  in  Figure  6,  we  picture  this  renewal  process 
with  renewals  occurring  at  replacement  every  t  units  of 


'  *•*  »*  ^  "'k  * 
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Fig.  6.  A  Block  Replacement  Renewal  Process 


time.  Thus,  in  this  model,  the  expected  renewal  length, 

E{L}  in  equation  (1.1),  equals  t.  To  obtain  V^,  the  cost 
t  h 

of  the  i  renewal  interval,  we  can  multiply  by  the 

number  of  failures,  N.,  in  the  i^^  renewal  interval 

i 

((i-l)t,  it]  and  then  add  in  order  to  account  for  the 
cost  of  the  single  replacement.  The  expected  cost  during 
a  typical  renewal  interval,  E{V},  can  be  expressed  as 

E{V}  =  M(t)  + 


where  M(t)  =E{N^},  the  renewal  function  [cf.  Ross  (1980: 
228-230)1 .  Thus,  assuming  the  replacement  process  will 
continue  indefinitely,  the  long  run,  average  cost  per 
unit  time,  B(t),  for  a  Block  Replacement  policy  is  a  func¬ 
tion  of  the  replacement  interval,  t,  and  can  be  expressed 
as 


B(t) 


E{V}  _ 
E{L} 


M(t)  +  C2 
t 


(1. 


As  an  aid  to  a  graphical  representation,  we  can  recast 
equation  (1.8)  in  terms  of  the  long  run,  average  stan¬ 
dardized  cost  per  unit  time,  i.e. 


M(t)  + 


B(t) 


=  B(t) 


(1.9] 


Solution  of  Block  Replacement 
When  M(t)  is  Specified 

If  M(t)  is  specified,  then  analytical  or  numerical 
solution  techniques  can  be  used  to  derive  the  optimal  t* 
which  minimizes  B(t).  However,  equation  (1.9)  has  a  use¬ 
ful  geometric  interpretation  when  we  plot  y  »  M(t)  versus 
t  as  depicted  in  Figure  7.  The  renewal  function,  M(t), 
determines  the  shape  of  this  plot  %diich  can  be  shown  to  be 
nondecreasing.  If  we  consider  the  curve  {t,  y»M(t)}  for 
0  £  t  £  oo,  the  slope,  b,  of  any  line  y  =  bt-C2/Cj^  from  the 
point  {0,-C2/Cj^}  to  {t,  M(t)}  equals  B(t).  The  tangent 
line,  as  drawn  in  Figure  7,  is  the  line  with  minimum 
slope  and  thus  represents  the  optimal  solution.  The  value 
t*  at  the  point  of  tangency  is  the  Block  Replacement  inter¬ 
val  which  achieves  a  minimum  for  B(t).  Therefore,  t*  is 
an  optimal  interval  with  B(t*)  equal  to  b*,  the  slope  of 
the  tangent  line.  The  utility  of  this  graphical  represen¬ 
tation  centers  on  a  visual  means  of  sensitivity  analysis 
for  various  values  of  With  a  straight  edge,  we  can 

easily  estimate  t*  for  different  cost  ratios. 

Depending  on  the  shape  of  the  renewal  function,  a 
finite  t*  may  not  exist.  In  the  case  of  exponential 
failures,  for  excunple,  the  renewal  function  is  linear,  and 


t* 


time,  t 


Fig .  7 .  Graphical  Representation  of  Block  Replacement 


thus  the  point  of  tangency  occurs  at  infinity  as  long  as 


>  0. 


Solution  to  Block  Replacement 
When  M(t)  is  Estimated 


In  Chapter  III,  we  estimate  the  optimal  Block 
Replacement  interval  using  the  above  graphical  technique 
and  an  empirical  estimate  of  the  renewal  function.  This 


empirical  renewal  function,  M^(t),  is  developed  from  the 


superposition  of  a  number  of  renewal  sample  paths. 
Suppose  we  have  identical  components  on  test  in  n  posi¬ 
tions,  and  we  replace  these  components  immediately  upon 


failure.  The  sequence  of  failure  times,  in  say  the 

.th 


i  cOTiponent  position,  represents  a  sample  path  of  the 
renewal  process.  As  in  Figure  8,  we  can  superimpose  the 
n  sample  paths  on  each  other  so  as  to  constitute  a  super¬ 


position  process  [cf .  Karlin  and  Taylor  (1975) ] .  The 

n 

superimposed  number  of  failures  E  N.(t),  form  the  series 

i=l  ^ 


{Z,,.,  Z,_,,  ...,  Z,  .}  where  the  last  superimposed 
(1)  (2)  (m) 
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Fig.  8.  Superposition  of  n  Renewal  Processes 

failure  time,  ^ (m) '  could  be  determined  by  some  life  test 
stopping  rule.  Then,  by  the  Strong  Law  of  Large  Numbers, 


^(1)  ^(2)  ^(3)  ^(4) 


lim 


£  N.(t) 


i*l 


n 


»  M(t)  . 


\  I 

There  are  two  methods  of  stopping  this  life  test, 
i.e.,  we  stop  at  a  specified  time,  t^,  or  stop  after  a 
certain  (say  r)  number  of  failures  have  occurred  at  each 
component  position.  In  the  first  case,  both  the  number 
of  superimposed  failures  and  their  failure  times  are 
random  variables.  Moreover,  the  number  of  superimposed 
failures  could  equal  zero.  In  the  latter  case,  only  the 
superimposed  failure  times  are  random  variables  because 
the  number  of  failures  is  fixed  at  m»nr.  In  the  inter¬ 
est  of  ensuring  some  superimposed  failure  times,  we  choose 


to  stop  after  r  failures  have  occurred  in  each  ccxnponent 
position. 


The  graphical  solution  technique  involves  esti¬ 
mating  t*  by  one  of  the  superimposed  failure  times,  z 
and  approximating  the  renewal  function  at  z^^j  as 


(j)' 


M(Z(.,) 


M  (z,  )  =  -1  . 

n'  (d)'  n 


Our  strategy  is  to  plot  y  =  Figure  7  and 

estimate  t*  from  the  point  of  tangency.  In  Chapter  III, 
we  investigate  (via  Monte  Carlo  simulation)  the  accuracy 
and  bias  of  our  estimator  as  well  as  the  relative  error  in 
the  objective  function  due  to  the  estimate.  We  conclude 
that  a  large  number  (greater  than  30)  of  component  posi¬ 
tions  may  be  necessary  in  order  to  reduce  the  li)celihood 
of  making  a  substantial  relative  error  in  B(t*)  due  to 
our  estimator.  If  the  underlying  distribution  is  exponen¬ 
tial  and  the  cost  ratio  small  (close  to  the 

origin) ,  there  is  a  large  probability  of  incorrectly 
estimating  a  finite  t*;  however,  the  relative  error  in 
B(t*)  will  tend  to  be  small.  We  propose  a  data  augmenta¬ 
tion  technique  for  use  with  scant  data  but  show  that  it 
may  adversely  affect  our  estimate  when  the  number  of  com¬ 
ponent  positions  is  large. 


Blind  Replacement  Model 
Kadncr  and  Jorgenson  (1962)  formulate  a  model 
where  a  replacement  occurs  at  intervals  of  time  t  +  K 


regardless  of  the  system's  state,  i.e.  if  we  are  uncertain 
about  the  system's  condition,  we  replace  at  periodic  time 
intervals  rather  than  inspect  and  replace  as  necessary. 
This  model  is  fundamentally  different  from  Age  and  Block 
Replacement  in  that  we  do  not  replace  at  failure  (an  impos¬ 
sible  action  if  we  do  not  know  the  instant  of  failure) . 
This  Blind  Replacement  strategy  is  reasonable  only  if  the 
cost  of  an  inspection  is  greater  than  the  cost  of  replace¬ 
ment  or  if  an  inspection  will  destroy  or  consume  the  item. 
Otherwise,  we  would  opt  for  a  Maintenance  Inspection 
strategy. 

Radner  and  Jorgenson  suppose  that  each  replacement 

takes  K  units  of  time  so  that  the  system  is  renewed  every 

t -f- K  units  of  time  as  shown  in  Figure  9.  Thus,  the 

expected  renewal  interval  length,  e{l}  from  equation  (1.1) 

equals  t -»■  K.  They  describe,  as  their  criterion,  "goodtime' 

(i.e.  availability) ,  G.  during  the  i^^  renewal  interval. 

;t 

Since  E{V}  =  ElG.)  =  J  R(x)dx,  the  long  run,  average 
availability  per  unit  time  is 

t 

f  R(x)dx 

rft)  -  E{V}  _  0*' 

"  Ein - t~K 


A  Blind  Replacement  Renewal  Process 


Solution  to  Blind  Replacement 


When  F(x)  is  Specified 

A  graphical  solution  to  Blind  Replacement  can 

be  found  by  plotting  the  curve  {t,  y=  J  R(x)dx}  for 

0 

0  £  t  £  00  as  in  Figure  10.  The  slope,  g,  of  any  line 

y  =  g(t+K)  from  the  point  {-K,0}  to  {t,  J  R(x)dx}  equals 

0 

the  value  of  the  objective  function  G(t)  at  t.  The 
tangent  line,  as  drawn  in  Figure  10,  is  the  line  with 
maximum  slope  and  thus  represents  the  optimum  solution. 

The  value  t*  at  the  point  of  tangency  is  the  Blind  Replace¬ 
ment  interval  which  achieves  a  maximum  for  G(t).  There¬ 
fore,  t*  is  an  optimal  interval  with  G(t*)  equal  to  g*, 
the  slope  of  the  tangent  line. 


Fig.  10.  Graphical  Representation  of  Blind  Replacement 

This  graphical  approach  is  similar  to  Bergman's 

method  shown  earlier  in  Figure  3.  However,  there  are 

several  minor  differences.  First,  with  this  approach, 

it  is  not  necessary  to  scale  the  conditional  mean  life, 
t 

f  R(x)dx).  Also,  the  value  of  G(t*)  can  be  obtained 


directly  from  the  picture  as  the  slope,  g*.  Finally, 
the  optimal  Blind  Replacement  interval  t*  can  be  read 
directly  on  the  abscissa. 


Solution  to  Blind  Replacement 
When  F(t)  is  Estimated 

We  suggest  estimating  the  optimal  Blind  Replace** 
ment  interval  t*  by  one  of  the  observed  failure  times  and 
approximating  the  conditional  mean  life  by  the  total  time 
on  test  statistic  divided  by  n.  In  Chapter  IV,  we  report 
on  the  results  of  Monte  Carlo  simulations  of  this  tech¬ 
nique  for  several  OFR  and  IFR  Weibull  distributions. 

We  conclude  that  the  total  time  on  test  statistic  produces 
positively  biased  estimates  of  t*  for  a  small  to  moderate 
replacement  time,  K.  If  we  use  the  survivor  time  on  test 
statistic  in  lieu  of  the  total  time  on  test  statistic, 
then  the  technique  produces  negatively  biased  estimates  of 
t*  for  moderate  to  large  K.  As  with  Age  Replacement,  we 
recommend  a  composite  estimator  for  Blind  Replacement  in 
order  to  mitigate  bias. 

Confidence  Intervals  for 
Maintenance  Planning 

In  the  graphical  solutions  to  the  three  maintenance 
strategies  discussed,  we  have  only  point  estimates  of  the 
optimal  replacement  interval  t*.  We  show,  in  Chapters  II, 
III,  and  IV,  that  these  graphical  techniques  are  biased 
and  that  the  degree  of  bias  diminishes  as  the  S€unple  size 
increases.  A  useful  measure,  in  addition  to  a  biased  point 


estimate  of  t*,  is  a  (1-a)  100  percent  confidence  interval 
for  t*.  In  Chapter  V,  we  propose  two  confidence  bounds 
techniques  based  on  Efron's  (1982)  bootstrap  resampling 
approach.  We  show,  via  Monte  Carlo  simulation  for  Age 
and  Blind  Replacement,  that  these  techniques  can  produce, 
in  many  cases,  useful  confidence  bounds.  However,  for 
Age  Replacement,  the  confidence  bounds  may  not  contain 
t*  when  there  are  very  large  or  very  small  differences 
between  the  cost  of  failure  and  the  cost  of  replacement. 

Summary 

Statistical  inference  issues  are  of  concern  when¬ 
ever  samples  are  used  to  infer  the  characteristics  of  a 
population.  Since  graphical  solution  techniques  for  the 
above-mentioned  maintenance  models  employ  estimates  of  the 
survival  probability  and  the  conditional  mean  life  based 
on  Scimples,  the  statistical  characteristics  of  these  esti¬ 
mates  are  of  interest. 

We  investigate  the  statistical  inference  issues  of 
graphical  solutions  for  Age  Replacement,  Block  Replacement 
and  Blind  Replacement  in  Chapters  II,  III,  and  IV  respec¬ 
tively.  In  many  (but  not  all)  cases,  we  observe  a  posi¬ 
tive  bias  in  the  estimators  for  t*  for  all  three  models. 
This  positive  bias  in  Age  and  Blind  Replacement  can  be 
offset  in  some  cases  with  a  composite  estimator  based  on 
the  total  time  on  test  statistic  and  the  survivor  time  on 


test  statistic. 
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Since  these  graphical  techniques  provide  only 
point  estimates  of  the  optimal  replacement  interval,  we 
are  also  interested  in  a  confidence  interval  methodology. 
In  Chapter  V,  we  propose  such  a  methodology  using  a  boot- 
strap  resampling  approach  and  demonstrate  its  utility  for 
both  the  Age  Replacement  and  Blind  Replacement  models. 


.AraC.-.W. 
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CHAPTER  II 


AGE  REPLACEMENT 

In  an  Age  Replacement  maintenance  scenario,  a 
single-part  item  is  replaced  or  renewed  to  good-as-new 
condition  at  failure  or  at  age  t  whichever  event  occurs 
first.  Barlow  and  Proschan  (1965)  provide  a  comprehensive 
synopsis  of  the  literature  relating  to  Age  Replacement  and 
a  concise  survey  of  useful  renewal  theory.  Their  develop¬ 
ment  of  a  transcendental  solution  for  the  optimal  Age 
Replacement  interval,  t*,  presupposes  that  the  costs  of 
failure  (C^^)  and  replacement  (C2)  are  known  with 
that  the  survival  distribution  belongs  to  the  Increasing 
Failure  feate  (IFR)  class,  and  that  any  failure  (or  life 
until  age  t)  is  instantly  replaced  the  moment  it  occurs. 
This  last  assumption  requires  instantaneous  knowledge  of 
condition  and  age,  information  which  may  be  very  expensive 
to  obtain.  This  expense  is  not  considered  in  their  formu¬ 
lation. 

Assuming  F(t)  is  differentiable  in  equation  (1.2), 
Barlow  and  Proschan  show  that  the  t*  which  minimizes  long 
run,  average  cost  per  unit  time  must  satisfy 


where  f(t)  is  the  derivative  of  F(t).  So  long  as 

F(t)  is  IFR,  they  note  that  equation  (2.1)  implies 
00 

t*  >  {C_/C, }  f  R(x)dx,  i.e.  we  would  never  want  to 
0 

schedule  a  preventive  replacement  prior  to  this  fraction 
of  the  mean  life  (which  is  not  the  same  as  the  conditional 
life) •  They  also  remark  that  t*  will  be  finite  if  the 
ratio  of  the  standard  deviation  to  the  mean  of  the  sur¬ 
vival  distribution  is  less  than  {1- (C2/Cj^)  } . 

If,  on  the  other  hand,  the  survival  distribution 
belongs  to  the  Decreasing  Failure  Rate  (DFR)  class  where 
the  instantaneous  probability  of  failure  decreases  with 
age,  then  Barlow  and  Proschan  point  out  that  Age  Replace¬ 
ment  is  inappropriate  since  a  used  item  has  greater 
expected  remaining  life  than  a  new  item  and,  hence,  there* 
is  no  incentive  to  replace.  As  a  special  case,  the  expo¬ 
nential  distribution  is  both  IFR  and  DFR;  its  failure  rate 
neither  increases  nor  decreases  but  rather  is  constant. 

We  can  easily  show  that  equation  (2.1)  has  no  solution  if 
the  survival  distribution  is  exponential.  Therefore,  we 
would  elect  to  forego  preventive  maintenance  in  this  case 
and  establish  the  optimal  Age  Replacement  interval  at 
infinity. 

The  IFR  class  contains  many  different  distribu¬ 
tions  including  the  truncated  normal  and  most  of  the  many 
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Weibull  and  gamma  distributions.  For  these  three  types 
of  distributions,  Glasser  (1967)  provides  mnemonic  graphs 
based  on  equation  (2.1)  to  aid  the  practitioner  in  deter¬ 
mining  the  optimal  Age  Replacement  interval.  Of  course, 
an  analyst  could  solve  equation  (2.1)  via  some  search 
technique,  or  for  that  matter,  search  techniques  could  be 
used  directly  on  the  objective  function,  in  equation  (1.2), 
should  F(t)  not  be  differentiable.  Nonetheless,  the 
failure  distribution  must  be  specified  a  priori.  In  a 
situation  with  scant  failure  data,  however,  there  may  be 
some  uncertainty  as  to  the  appropriate  failure  distribu¬ 
tion.  With  this  problem  in  mind,  Arunkumar  (1972)  con¬ 
siders  a  nonparametric  approach  to  estimating  F(t). 

Arunkumar  (1972)  introduces  the  concept  of  an 
empiric  approximation  to  the  survival  distribution  in  Age 
Replacement  without  recommending  a  plotting  position.  He 
shows  that  an  estimator  of  the  optimal  Age  Replacement 
interval,  based  on  some  empirical  survival  distribution, 
is  strongly  consistent.  In  an  applied  innovation,  Bergman 
(1977a)  expands  on  this  idea.  He  uses  an  empirical  cdf 
based  on  the  California  plotting  position  as  an  approxi¬ 
mation  to  the  failure  distribution.  As  we  discuss  in 
Chapter  I,  this  approximation  leads  to  the  total  time  on 
test  statistic  which,  when  divided  by  n,  approximates  the 
conditional  mean  life.  In  this  chapter,  we  consider  the 
rapidity  by  which  Bergman's  method,  which  we  outlined  in 


Chapter  I,  converges  to  the  optimal  Age  Replacement  inter¬ 
val. 

Although  we  focus  on  Barlow  and  Proschan's  Age 
Replacement  model,  there  are  variations  on  this  model  which 
also  may  be  of  interest  to  the  maintenance  community  in 
terms  of  graphical  solutions.  We  do  not  pursue  this  issue 
here,  but  rather,  we  detail  several  variations  to  indicate 
a  potential  for  further  research.  For  example,  Scheaffer 
(1971)  supposes  that  the  cost  function,  C(t),  includes  a 
term  which  increases  with  an  item's  age  (to  account  for 
increasing  replacement  cost  "due  to  depreciation  or  wear"). 
His  results  require  explicit  knowledge  of  the  relationship 
between  cost  and  age,  which  may  be  exceedingly  difficult 
to  obtain  in  practice.  Weiss  (1956)  formulates  an  Age 
Replacement  model  which  includes  the  time  and  cost  of 
repair  as  well  as  the  cost  of  failure  and  replacement.  He 
indicates  that  an  optimal  solution  may  not  always  exist. 


Graphical  Solutions  for  Age  Replacement 


The  total  time  on  test  statistic  is  used  by 
Bergman  (1977a)  to  develop  a  graphical  solution  for  Age 
Replacement.  The  total  time  on  test  statistic. 


when  divided  by  the  total  time  on  test,  T(tj^j),  is  the 
scaled  total  time  on  test,  U^,  at  age  The  plot  of 
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versus  i/n  is  the  total  time  on  test  (TTT)  plot. 

Barlow  and  Campo  (1975)  point  out  that,  as  n-^co,  an  IPR 
distribution  results  in  a  concave  TTT  plot  on  the  inter¬ 
val  (0,1),  whereas  a  distribution  from  the  DPR  class  will 
have  a  convex  plot.  An  exponential  life  distribution 
(which  belongs  to  both  the  IPR  and  DPR  classes)  will  have 
a  linear  TTT  plot  passing  through  the  origin  at  a  45 
degree  angle,  i.e.  the  plot  is  both  concave  and  convex. 

An  empirical  life  distribution  need  not  belong  to  either 
the  IPR  or  DPR  class,  and  thus,  its  TTT  plot  would  not 
necessarily  be  concave  or  convex.  Barlow  (1979)  further 
describes  properties  of  TTT  plots  and  shows  that  concave 
(convex)  TTT  plots  imply  IPR  (DPR)  survival  distributions. 
If  the  survival  distribution  belongs  to  the  Increasing 
Failure  Rate  Average  (IFRA)  class  or  to  the  Decreasing 
Failure  Rate  Average  (DFRA)  class,  then  he  proves  that 
the  TTT  plot  will  be  anti-starshaped  for  IFRA  distribu¬ 
tions  and  starshaped  for  DFRA  distributions.  A  starshaped 
TTT  plot  is  one  that  is  on  or  below  a  linear  TTT  plot 
passing  through  the  origin  at  a  45  degree  angle  but  that 
is  not  necessarily  convex.  Conversely,  an  anti-starshaped 
TTT  plot  is  on  or  above  a  45  degree  line  through  the 
origin  but  is  not  necessarily  concave.  Although  an  anti¬ 
starshaped  TTT  plot  is  not  sufficient  evidence  that  the 
life  distribution  is  IFRA,  Barlow  (1979)  remarlts  that  this 
is  sufficient  to  conclude  that  a  survival  distribution 
"is  at  least  NBUE  [New  Bettor  than  Used  in  Expectation, 
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a  class  of  distributions  which  includes  the  IFRA  class] 
if  not  IFRA." 

Berman's  method  employs  the  TTT  plot  as  a  scaled 
representation  of  the  conditional  mean  life,  i.e.,  is 
proportional  to  T(t^j^j)/n.  Hence,  the  ^ which  maxi¬ 
mizes  the  ratio  U^/ { (C2/ ^  ^  will  also  maxi¬ 

mize  the  ratio  (T  (t  ) /n)  /  { (Cj/ (€^^-€2) )  +  (i/n)  } .  As 
Bergman  (1977a)  and  others  note,  if  tjjj  “^(n) 
seems  reasonable  to  estimate  the  age  replacement  interval 
as  infinity,"  i.e.  t*  =  «».  This  rationale  is  consistent 
with  the  Age  Replacement  solution,  t*  =  «>,  should  the  under¬ 
lying  survival  distribution  be  exponential  or  DFR.  If 
the  survival  distribution  is  exponential,  then  the  TTT 
plot  should  approximate  a  45  degree  line  from  the  origin. 

If  it  is  DFR,  then  the  TTT  plot  should  approximate  a  convex 
function.  In  either  case,  the  point  of  tangency  would 
likely  correspond  with  provided  that  there  was  ample 

data  to  give  a  good  approximation. 

With  small  sample  sizes,  however,  it  may  be  diffi¬ 
cult  to  conclude  that  the  data  was  generated  by  an  expo¬ 
nential  or  DFR  distribution.  Given  some  failure  data,  we 
may  be  able  to  apply  a  variety  of  goodness-of-f it  tests 
to  reject  the  null  hypothesis  that  our  data  could  have 
come  from  an  exponential  (or  DFR)  distribution.  However, 
if  we  can  not  reject  the  null  hypothesis,  then  we  still 
can  not  confirm  that  the  data  are  indeed  exponential  (or 
DFR)  since  it  could  be  that  our  sample  size  is  insufficient 


.  j.  .r,  «•, 


to  detect  that  the  null  hypothesis  is  false.  As  an 
alternative  goodness-of-f it  test,  Barlow  and  Campo  (1975) 
propose  counting  the  number  of  times  the  empiric  TTT  plot 
crosses  the  exponential  TTT  plot.  Presumably,  exponen¬ 
tial  observations  will  result  in  an  empirical  TTT  plot 
that  crosses  the  45  degree  line  from  the  origin  a  large 
number  of  times.  Bergman  (1977b)  provides  the  exact  and 
asymptotic  distributions  for  the  number  of  crossings  for 
exponential  data.  When  using  Bergman’s  graphical  method 
for  Age  Replacement  without  knowledge  of  the  underlying 
distribution,  one  should  be  concerned  about  the  proba¬ 
bility  of  estimating  a  finite  Age  Replacement  interval 
when,  in  fact,  the  distribution  is  exponential  (or  DFR) . 

To  address  this  issue,  we  simulated  5000  applica¬ 
tions  of  Bergman's  method  to  IFR  and  DFR  lifetime  data  for 
various  sample  sizes.  As  one  would  expect,  the  expected 
value  of  the  longest  lifetime,  E{T^^j},  increases  for 
larger  sample  sizes  from  exponential  and  DFR  distributions 
As  we  show  in  Table  1,  Bergman’s  technique  when  applied  on 
exponential  observations  resulted  in  about  a  50  percent 
probability  of  incorrectly  estimating  a  finite  Age  Replace 
ment  interval  when  the  cost  ratio  C^/  (Cj^-C2)  equals  1.0. 

In  this  table,  is  the  probability  of  making  the 

correct  decision  not  to  have  preventive  replacement.  This 
probability  is  sensitive  to  the  cost  relationship.  If 
^2^  is  greater  than  one  (suggesting  only  a  nominal 

difference  between  the  cost  of  failure  and  the  cost  of 


replacement) ,  then  there  is  a  higher  probability  of  reach¬ 
ing  the  correct  conclusion  of  no  preventive  maintenance. 

On  the  other  hand,  if  Q^t  is  close  to  the  origin 

(as  would  be  the  case  with  a  substantial  difference  in 
costs) ,  then  there  is  a  large  probability  of  reaching  the 
incorrect  conclusion  of  a  finite  Age  Replacement  interval. 
The  reason  for  this  phenomenon  is  that  close  to  the  origin 
there  is  a  "seunple  clutter"  effect.  This  occurs  when  a 
tangent  line  (passing  through  [-C^/ (C^-C2  »0]  )  strilces  a 
point  on  the  TTT  plot  associated  with  a  failure  time, 
t^jj,  less  than  where  the  tangent  point  may  be  only 

slightly  above  the  45  degree  line.  Thus,  for  very  large 
differences  between  and  C2/  the  practitioner  ought  to 
view  a  finite  estimate  of  t*  with  caution  regardless  of 
the  sample  size.  A  prudent  course  of  action  would  be  to 
checlt  for  exponentiality  (or  DFR)  with  some  statistical 
test  and  accept  the  finite  estimate  only  if  the  data  are 
not  exponential  (or  DFR) . 

The  economic  consequences  of  ignoring  this  issue 
might  be  substantial.  If  the  observed  lifetimes  are 
exponential  with  scale  parameter  X,  then  t*  =  «>  and 
C('”)  =XCj^.  However,  if  Bergman's  technique  produces  a 
finite  estimate,  t*,  of  t*,  then  C(t*)>  C(t*).  In 
Table  2,  we  show  the  average  relative  error  in  the  objec¬ 
tive  luncLion,  (C ( t* ) -C ( t* ) ) /C ( t* ) ,  based  on  the  simula¬ 
tion  of  5000  replications  of  Bergman's  technique  on  expo¬ 
nential  observations.  Note  that,  as  the  cost  ratio 


OBJECTIVE,  C(t*),  DUE  TO  BERGMAN’S  TECHNIQUE 
FOR  OBSERVATIONS  FROM  AN  EXPONENTIAL 
DISTRIBUTION  (MEAN  »  200) 


decreases,  the  sample  clutter  effect  aggravates  the 

A 

average  relative  error  in  C(t*).  We  remark  that  this 
average  relative  error  is  within  a  nominal  10  percent 
range  for  the  range  of  cost  ratios  shown. 

For  the  Weibull  family  of  distributions,  the 
probability  of  estimating  an  infinite  replacement  age 
diminishes  for  increasingly  IFR  Weibull  distributions  and 
increases  for  more  DFR  Weibull  distributions  as  we  see 
in  Table  3.  Of  concern,  perhaps,  is  the  relatively  large 
probability  of  incorrectly  electing  to  forego  preventive 
replacement  (this  is  in  Table  3  for  items  with  a 

slightly  IFR  survival  distribution) .  In  other  words, 
with  Bergman's  method,  we  may  erroneously  call  for  the 
preventive  replacement  of  an  item  not  requiring  it  and, 
conversely,  in  no  preventive  replacement  for  an  item  sub¬ 
ject  to  wearout.  The  average  relative  error  in  C(t*)  for 
these  distributions  is  shown  in  Table  4. 

Bergman's  method  approximates  the  failure  dis¬ 
tribution,  F(t),  by  the  California  plotting  position,  i/n 
and  employs  the  total  time  on  test  statistic  to  approxi¬ 
mate  the  conditional  mean  life.  Moreover,  the  true, 
optimal  interval  is  approximated  by  one  of  the  observed 
failure  times.  The  performance  of  these  three  approxi¬ 
mations  is  of  interest  because,  as  we  show  below,  their 
synergistic  effect  on  the  estimator  can  be  biased.  We 
simulated  5000  replications  of  Bergman's  technique  on 
three  IFR  Weibull  distributions  with  different  shapes 
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which  result  in  a  finite  t*  for  any  cost  ratio, 

0  < -C^/ (C^-C^)  «»•  These  different  shaped  distributions 
are  convex  ordered,  and  hence,  a  large  shape  parameter 
distribution  is  said  to  be  more  IFR  them  a  small  shape 
parameter  distribution.  Accordingly,  these  three  dis¬ 
tributions,  which  we  consider,  are  nominally,  moderately, 
and  highly  IFR,  respectively.  For  increasingly  IFR 
distributions,  there  is  a  decreasing  probability  of  esti¬ 
mating  t*  as  infinity  (see  )  in  Table  5) .  We 

observe,  from  Table  5,  that  the  expected  value  of  Bergman's 
estimator,  E{T*},  has  a  positive  bias  for  moderate  to 
highly  IFR  Weibull  distributions  and  a  negative  bias  for 
the  nominally  IFR  distribution. 

As  an  alternative  to  estimating  t*  with  a  graphi¬ 
cal  technique  that  approximates  the  conditional  mean  life 
using  the  total  time  on  test  statistic,  we  could  replace 
this  statistic  with  the  survivor  time  on  test  statistic. 

The  resultant  estimate  has  a  negative  bias  for  the  low  to 
moderately  IFR  distributions  and  a  positive  bias  for  the 
highly  IFR  distribution  as  we  show  in  Table  6.  A  combina¬ 
tion  of  these  two  estimators,  then,  may  result  in  less 
bias  than  either  estimator  alone  if  the  failure  distribu¬ 
tion  is  moderately  IFR. 


Summary 

Graphical  solutions  to  the  Age  Replacement  problem 
offer  an  attractive  alternative  to  analytical  solutions 
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TABLE  6 


EXPECTED  VALUE  OF  BERGMAN'S  ESTIMATOR,  E{T*};  PROBABILITY 
THAT  t(j)  IS  THE  LAST  FAILURE  TIME,  P{t(n)};  AMD  AVERAGE 
%  RELATIVE  ERROR  (%  RE)  IN  C(t*)  FOR  OBSERVATIONS  FROM 
VARIOUS  WEIBULL  DISTRIBUTIONS  (MEAN  »  200) 

USING  STT  STATISTIC 


Sample 

Size 

for 

C2/(C^-C 

2)  -  1.0 

and  Weibull  Shape 

of 

1 

E{T*} 

.5 

%  RE 

2. 

E{T*} 

0 

%  RE 

3.0 

E{T*} 

%  RE 

10 

266.2 

6.52 

221.5 

3.35 

185.9 

2.77 

20 

309.6 

3.26 

235.1 

1.93 

185.6 

1.97 

30 

328.5 

2.35 

239.0 

1.49 

185.2 

1.60 

40 

342.2 

1.86 

242.4 

1.19 

185.2 

1.33 

50  / 

352.5 

1.55 

243.0 

1.04 

185.1 

1.15 

60 

363.2 

1.27 

244.4 

0.93 

184.3 

1.03 

70 

372.6 

1.11 

245.5 

0.86 

183.6 

0.91 

80 

376.9 

1.00 

245.2 

0.78 

183.7 

0.86 

90 

380.7 

0.93 

245.3 

1 

0.73 

183.5 

0.80 

100 

386.6 

0.81 

246.8 

0.70 

183.3 

0.74 

«0 

481.0 

- 

249.9 

- 

182.7 

- 

and  their  inherent  uncertainty  of  the  failure  distribution 
and  the  complexity  of  search  routines  to  solve  for  the 
minimizer  of  either  equations  (1.2)  or  (2.1).  Bergman 
(1977a)  has  proposed  a  simple  graphical  method  involving 
the  total  time  on  test  plot  and  an  empirical  failure  dis¬ 
tribution.  The  seumpling  issues  of  this  approach  are  of 
interest  when  we  consider  statistical  inferences  from  a 
set  of  observed  failures.  We  find  that  Bergman's  method 
may  incorrectly  lead  to  finite  estimates  of  the  optimal 
Age  Replacement  interval  if  there  are  substantial  differ¬ 
ences  between  the  costs  of  failure  and  replacement  and  if 
the  (unknown)  parent  failure  distribution  is  exponential 
or  DFR.  The  economic  consequences,  however,  in  terms  of 
the  percent  relative  error  are  nominal.  Nonetheless, 
Bergman's  method  may  accurately  estimate  optimal  Age 

Replacement  intervals  for  many  IFR  distributions  with 

/ 

samples  as  small  as  20  observations. 


BLOCK  REPLACEMENT 


Block  Replacement  is  a  maintenance  strategy  which 
does  not  require  knowledge  of  an  item's  age.  Such  infor¬ 
mation  may  be  expensive  or  impossible  to  obtain.  Rather, 
with  Block  Replacement,  we  replace  or  renew  to  good-as- 
new  condition  at  failure  and  at  periodic  intervals  of  time 
regardless  of  when  the  last  failure  occurred.  Barlow  and 
Proschan  (1965;  1975)  discuss  Block  Replacement  and  remark 
that,  when  compared  to  Age  Replacement,  Block  Replacement 
results  in  more  removals  regardless  of  the  survival  dis¬ 
tribution  and  in  fewer  failures  when  the  survival  dis¬ 
tribution  is  IFR.  Although  more  useful  life  is  discarded 
under  Block  Replacement  than  Age  Replacement,  there  is  no 
necessity  to  record  age. 

Presuming  a  failure  cost  and  a  preventive 
replacement  cost  C2>  Barlow  and  Proschan  (1965)  formu¬ 
late  a  long  run,  average  cost  per  unit  time  function, 

B(t),  as 


C,  M(t)  +  C- 


They  show  the  first  order  necessary  condition  to  be 


m(t)t  -  M(t) 


(3.21 


where  m(t)  is  the  derivative  of  the  renewal  function, 

M(t) ,  which  is  the  sum  of  the  k-fold  convolution  of  the 
failure  distribution  with  itself  for  k>l,  As  we 

explain  in  Chapter  I,  the  renewal  function  gives  the 
expected  number  of  failures  in  the  interval  (0,t] .  Barlow 
and  Proschan  reason  that  B(t)  approaches  infinity  as  t 
approaches  zero  and  conclude  that  B(t)  has  a  minimizer, 
t*,  in  the  interval  (0,<»]  where  t*<»  e»  is  interpreted  as 
replacement  at  failure  only.  As  we  also  explain  in 
Chapter  1,  a  graphical  clarity  can  be  obtained  by  consider' 
ing  the  standardized  objective 


B(t)  = 


B(t) 

^1 


M(t)  +  ^ 
^1 


(3.3 


The  renewal  function  can  be  expressed  in  terms  of 
the  fundamental  renewal  equation 


M(t)  =  F(t)  + 


/ 


t 

M(t-x)dF(x) 


(3.4 


which  is  often  solved  via  the  Laplace-Stieltjes  transform. 
There  are  several  basic  results  for  selected  processes 
which  we  use  in  later  examples.  Karlin  and  Taylor  (1975) 


cheuracterize  a  Poisson  process  with  pareuneter  X  as  a 
renewal  process  and  note  that  M(t)  *  Xt,  a  linear  func¬ 
tion.  They  also  give  a  closed  form  solution  for  the 
renewal  function  of  an  IFR  gamma  distribution  with  shape 
parameter  of  two  and  scale  parameter  beta  as 


This  renewal  function  is  a  convex  function  in  t.  As  for 
the  DFR  class  of  distributions.  Brown  (1980)  proves  that 
their  renewal  functions  are  concave.  We  make  use  of  these 
geometric  properties  in  a  graphical  solution  to  Block 
Replacement . 

Earlier  in  Figure  7,  we  depicted  a  geometric  solu¬ 
tion  to  the  Block  Replacement  model  by  plotting  y^MCt) 
versus  t.  If  the  lifetime  distribution  is  exponential, 
then  the  renewal  function,  M(t) ,  is  a  linear  function  in 
t  so  that  there  is  no  finite  temgent  point  as  long  as 
C2/Ci>0.  Hence,  the  optimal  Block  Replacement  policy 
for  exponential  lifetimes  is  to  replace  at  failure  only. 
However,  an  IFR  gamma  distribution,  with  shape  parameter 
two  and  scale  parameter  beta,  has  a  convex  renewal  func¬ 
tion  in  t.  Barlow  and  Proschan  show  that  this  renewal 
function  results  in  a  finite  optimal  Block  Replacement 
interval  so  long  as  and  C2  satisfy  the  inequality 
C2/C1 <  .25.  DFR  distributions,  on  the  other  hand,  have 
concave  renewal  functions  and,  hence,  optimal  Block 


Replacement  intervals  of  infinity.  Notwithstanding  these 
geometric  properties,  an  analytic  or  numerical  solution  to 
Block  Replacement  requires  specification  of  the  renewal 
function  or  at  least  the  lifetime  distribution.  As  we 
outlined  in  Chapter  I,  we  can  address  our  uncertainty  of 
the  renewal  function  by  approximating  it  by  an  empirical 
renewal  function  developed  from  a  superposition  process 
as  in  Figure  8. 

n 

The  superimposed  number  of  failures,  £  N.(t) 

i*l  ^ 

form  the  series  ^  ^2) '  *•*'  ^  (m)  ^  where  the  last 

superimposed  failure  time,  ^  determined  by  some 

arbitrary  life  test  termination  time  t,  i.e.  Z £  t. 
Then,  by  the  law  of  large  numbers. 


lim  i=l 


=  M(t) 


For  n  finite,  we  can  approximate  M(t)  as 


M(t)  ::  M^(t)  =  ^  . 


The  number  of  events,  m,  in  this  superposition  process  is 
a  function  of  t.  There  are  two  methods  of  determining  t 
which  are  akin  to  the  notion  of  life  test  censoring. 

Either  the  life  test  is  stopped  at  a  specified  time,  t^, 
or  it  is  terminated  after  a  particular  number  of  failures. 


say  "r,"  have  occurred  at  each  component  position.  In  the 
first  case,  both  the  number  of  superimposed  failures 
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(possibly  zero)  and  theijE'  failtil^e  times  are  random  vari¬ 
ables.  In  the  latter  case,  only  the  super in^osed  failure 
times  are  random  variables  since  the  number  of  failures  is 
fixed  (at  some  positive  quantity) .  We  choose  to  terminate 
after  r  failures  occur  in  each  consonant  position,  thereby 
guaranteeing  a  quantity  of  m  superinqposed  failures.  Here, 
m  will  be  equal  to  the  product  of  the  number  of  component 
positions,  n,  eind  the  nximber  of  failures  at  each  position, 
r . 

The  order  in  which  these  r  lifetimes  will  occur  in 
any  component  position  is  r2mdom.  For  any  given  set  of  r 
lifetimes,  there  are  r!  possible  sample  paths  representing 
the  permutations  in  which  these  lifetimes  could  occur. 

We  can  make  full  use  of  the  information  in  the  n  sets  of 
r  lifetimes  by  sampling  within  each  set  to  obtain  more 

sets  of  r  lifetimes  with  a  different  order  of  occurrence. 

/ 

We  discuss  this  issue  later  as  a  data  augmentation  tech¬ 
nique  . 

.1 

Graphical  Solution  to  Block  Replacement 

Following  an  argument  similar  to  that  of  the 

graphical  approach  to  Age  Replacement,  we  can  solve  an 

approximation  to  equation  (3.3)^  by  direct  numeration 

methods.  Given  a  set  of  super inqposed  failure  tiroes, 

% 

{z^j,  ...,  we  estimate  the  roinimizer,  t*,  of 

B(t)  by  one  of  the  summed  failure  times,  t*sz^jj,  such 
that  .  ^ 
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B(t*)  :  = 


(j) 


(3.6) 


Since  = j/n,  we  can  further  approximate  B(Zjjj) 

by  this  empiric  renewal  function  such  that 


B(z 


(3.7) 


Our  strategy  is  to  Initiate  a  life  test  with  n  conponent 
positions  and  run  the  test  until  r  failures  have  occurred 
in  each  component  position.  We  then  evaluate  equation 
(3.7)  for  each  summed  failure  time,  sind  identify 

that  z^jj  =  t*  which  minimizes  The  accuracy  of 

our  estimate  will,  in  part,  be  a  function  of  how  many 
component  positions  we  use  in  the  life  test  as  well  as 
the  number  of  failures  at  each  position.  If  t*  happens  to 
be  then  two  possibilities  arise;  either  the  true 

optimal  Block  Replacement  interval  is  indeed  infinity  or 
the  optimal  interval,  t*,  is  greater  than  z  but  still 
finite.  To  resolve  this  question,  a  practitioner  should 
extend  the  termination  point  (i.e.  increase  r)  by  some 
reasonable  amount  with  more  life  testing  in  each  component 
position. 

We  investigated  the  question  of  accuracy  by  simu¬ 
lating  1000  applications  of  the  above  methodology  for  IFR 
and  DFR  lifetime  data  with  various  numbers  of  component 


positions.  In  Table  7,  we  show  the  expected  value  of  our 
estimator,  E(T*},  for  increasing  numbers  of  component 
positions,  n,  when  the  survival  distribution  is  an  IFR 
gamma  distribution  (shape  parameter  two,  scale  parameter 
one) .  Recall  that  for  cost  ratios  (€2/0^)  less  than  .25, 
there  should  be  a  finite  Block  Replacement  interval.  As 
we  see  in  Table  7,  the  expected  value  of  our  estimator, 
E{T*},  overestimates  t*  for  large  differences  between 
and  C2,  say  ^  converges  toward  t*  as  the 

number  of  positions  increase.  We  note  that  the  proba¬ 
bility  of  t*  being  the  last  observed  superimposed  failure, 

P{t*  =  z,  .}  (which  we  denote  as  P{z  })  ,  is  small  as  long 
(m)  m 

as  the  cost  ratio  is  small.  However,  as  the  cost  ratio 
increases  there  is  a  greater  likelihood  that  t*  «  z  . 
Thus,  for  cost  ratios  as  small  as  ten  to  one,  (which 

A 

should  result  in  a  finite  t*  for  this  distribution)  there 

/ 

is  a  substantial  probability  of  incorrectly  estimating  the 
optimal  Block  Replacement  interval  as  infinity.  This 

A 

accounts  for  the  lack  of  convergence  of  E{T*},  as  shown 
in  Table  7,  when  t*  «  “  (as  would  be  the 

case  if  the  point  of  tangency  corresponds  to  then 

by  the  elementary  renewal  theorem,  B(t*  =  '»)  =Cj^/y  where  y 
is  the  expected  value  of  a  lifetime  generated  by  the 

A 

underlying  distribution.  For  finite  t*,  we  can  compute 
B(t*)  provided  that  we  have  a  closed  form  solution  for 

A 

M(t*)  as  in  equation  (3.5).  Accordingly,  in  Table  8,  we 
show  the  average  relative  error,  { [B (t*) -B (t*) ] /B ( t*) } , 


in  the  objective  function  due  to  our  estimator  T*  when  the 
failure  distribution  is  gamma  with  shape  parameter  equal 
two  and  scale  equal  one.  These  results  lead  to  the  intui¬ 
tive  conclusion  that  when  there  is  a  large  cost  differ¬ 
ence  between  and  C21  we  should  employ  this  solution 
technique  with  a  large  number  of  component  positions. 

Renewal  function  plots  of  DFR  distributions  are 
concave.  The  exponential  distribution,  as  a  special  case 
DFR  distribution,  has  a  renewal  function  that  is  linear 
in  t.  Theoretically,  as  long  as  the  cost  ratio  ^2/0^  is 
finite,  the  slope  b*  of  the  tangent  y  =  b*t-C2/Cj^  (shown 
in  Figure  7)  will  be  zero  and  t*  =  <».  However,  with  an 
empirical  renewal  function,  we  are  approximating  this 
linear  (or  concave)  renewal  function  with  a  piece-wise 
linear  function  which  may  have  a  finite  point  of  tangency. 
In  fact,  as  we  show  in  Table  9,  smaller  cost  ratios  (i.e. 

t 

closer  to  the  origin) ,  result  in  decreasing  probability  of 
correctly  estimating  an  infinite  Block  Replacement  inter¬ 
val.  However,  the  consequences  of  this  statistical  error 
appear  to  be  minimal  as  evidenced  by  the  low  relative 
error  in  B(t*)  shown  in  Table  9. 

Data  Augmentation  Technique 
for  Block  Replacement 

The  empirical,  renewal  function  estimator  dis¬ 


cussed  above  may  require  failure  data  from  a  leurge  number 
of  component  positions  to  estimate  adequately  the  optimal 
Block  Replacement  interval,  t*,  of  an  IFR  class 


distribution.  For  small  numbers  of  component  positions, 
the  estimator  may  overestimate  t*.  The  impact  of  scant 
data  might  be  diminished,  however,  by  "creating"  more  data 
with  a  sampling  scheme.  In  this  scheme,  the  r  lifetimes 
in  a  percentage  of  the  component  positions  are  randomly 
reordered  and  treated  as  additional  component  position 
failure  data.  These  additional  data,  of  course,  are  not 
independent  of  the  original  failure  data.  However,  given 
a  large  number  of  lifetimes  in  each  position,  a  randomly 
selected  sequence  of  failure  times  may  result  in  another 
saunple  path  which  is  weakly  dependent  on  the  original 
sample  path.  Hence,  our  estimating  technique  might  be 
robust  insofar  as  the  requirement  for  independent  sample 
paths  is  concerned. 

To  illustrate  this  seimpling  idea  via  Monte  Carlo 
simulation,  we  used  a  sequential  sampling  plan  where  we 
randomly  sampled  the  r  lifetimes  in  alternate  component 
positions.  We  drew  r  observations  with  replacement  and 
recorded  them,  in  the  order  drawn,  as  failures  in  pseudo¬ 
positions.  Then,  we  superimposed  the  failure  times  in 
both  the  actual  positions  and  the  pseudo-positions  to  form 
the  series  ty  £  •  •  •  £  Y  }  where  m  now  equals  1.5  nr. 
With  this  seunpling  plan,  we  artificially  increased  the 
number  of  superimposed  data  points  by  50  percent  and  (we 
hope)  improved  our  optimal  Block  Replacement  estimator. 

In  Table  10,  we  display  the  results  of  1000  iterations  of 
this  data  augmentation  technique  to  produce  an  expected 
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value  of  the  estimate  t'  compared  to  the  expected 

/S 

value  of  the  previous  estimate  t*  for  a  gamma  (shape  2, 
scale  1)  distribution  with  various  cost  ratios  resulting 
in  a  finite  t*. 

At  first  glance,  this  data  augmentation  technique 
does  not  appear  to  result  in  an  improved  estimator  except 
when  the  cost  ratio  is  very  small  (see 

Table  10)  and  there  are  few  numbers  of  component  positions. 

A. 

However,  when  the  relative  error  in  B{t*}  is  compared  to 

A 

the  relative  error  in  B{t*},  shown  in  Table  10,  we  note 

A 

that  the  estimate  t'  results,  on  average,  in  smaller  rela¬ 
tive  error.  This  advantage  diminishes  as  the  number  of 
component  positions,  n,  increase,  and  for  n  greater  than 
30  (in  Table  10) ,  t*  would  be  preferred  to  t’ .  Notwith¬ 
standing  this  disadvantage,  our  data  augmentation  tech¬ 
nique  may  be  useful  when  there  are  scant  data. 


Summary 

A  graphical  solution  to  Block  Replacement  can  be 
developed  by  using  an  empiric  estimator  of  the  renewal 
function  based  on  a  superposition  process  concept.  The 
idea  is  to  superimpose  r  lifetimes  from  each  of  n  ccan- 
ponent  positions  as  ^  1  1  ^  }  and  estimate  the 

expected  number  of  renewals  in  the  interval  [0,z^^j]  as 
i/n.  The  graphical  technique,  then,  identifies  that 
which  corresponds  to  the  point  of  tangency  on  a  plot  of  the 
empiric,  renewal  function  for  a  line  drawn  through  the 
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point  tO,-C2/Cj^].  This  is  an  estimator  of  the 

optimal  Block  Replacement  interval,  t*,  which  minimizes 
long  run,  average  cost  per  unit  time.  The  accuracy  of 
this  estimator  can  be  enhanced  by  increasing  the  number 
of  ccxnponent  positions.  This  is  especially  important 
when  the  difference  between  the  cost  of  replacement  and 
the  cost  of  failure  is  great. 


CHAPTER  IV 

BLIND  REPLACEMENT 

"Inspection  policies"  are  a  class  of  maintenance 
models  characterized  by  two  kinds  of  uncertainty.  For 
this  class,  not  only  are  we  uncertain  of  an  item's  life¬ 
time,  which  can  be  expressed  in  terms  of  a  probability 
distribution,  but  we  are  also  uncertain  of  the  item's 
condition  (good  or  failed).  Savage  (1956)  and  others 
have  referred  to  this  class  as  that  of  "preparedness 
models"  pertaining  to  items  in  storage  that  may  be  called 
into  service  during  some  random  emergency.  We  choose  to 
broaden  this  definition  to  include  in-service  items,  with 
an  emphasis  on  the  matter  of  item  condition  and  its  uncer' 
tainty.  To  determine  condition,  we  generally  inspect  and 
then  replace  (or  repair)  a  defective  item  upon  discovery. 
For  some  items,  the  cost  of  inspection  exceeds  the  cost 
of  replacement.  For  example,  automotive  engine  crankcase 
lubricant  is  commonly  replaced  at  periodic  intervals 
rather  than  inspected  for  condition  and  replaced  as  neces 
sary.  In  this  case,  we  presume  the  cost  of  inspection  (a 
laboratory  analysis)  greatly  exceeds  the  cost  of  replace¬ 
ment.  For  other  items  such  as  explosive  devices,  an 
inspection  involves  consumption  of  the  item,  and  hence 


the  mean  exponential  lifetime  increases.  Regardless  of 
the  lifetime  distribution,  we  note  that  to  use  either 
equation  (4.1)  or  (4.2),  we  must  specify  the  survival  dis¬ 
tribution,  R(x) ,  as  well  as  the  replacement  time,  K.  in 
practice,  we  may  be  uncertain  of  either  and,  therefore, 
may  find  empirical  techniques  similar  to  those  of  Bergman 
to  be  useful. 

Graphical  Solutions  for  Blind  Replacement 
There  are  two  straightforward  geometric  solution 
techniques  for  Blind  Replacement  which  use  the  same 
approximations  that  Bergman  uses  in  Age  Replacement  and 
which,  as  a  consequence,  may  be  subject  to  bias.  Given  a 
set  of  complete,  life  test  data,  (i)  ^  i  ^  (n)  ^ 
can  evaluate  either  equation  (4.1)  or  (4.2)  for  each 
(i  =  l,  ...,  n)  and  estimate  t*  by  the  maximizer, 
from  the-'  set  of  failure  times.  Using  equation  (4.1),  we 
can  numerically  evaluate  the  objective  function  at  each 
t^j^j  and  identify  the  which  maximizes  G(t^j^j)  as  an 

estimate,  t*,  of  the  optimal  interval.  We  refer  to  this 
as  a  "direct  numeration"  method.  Another  method  evaluates 
equation  (4.2)  for  each  and  it  identifies  that 

which  most  closely  satisfies  the  quality  as  the  "first 

fl 

order"  estimate,  t  ,  of  the  optimal  Blind  Replacement 


interval • 


Direct  Numeration  Solutions 

The  graphical  argument  for  the  direct  numeration 

solution  was  portrayed  earlier  in  Figure  10.  Our  approach 

involves  estimating  t*  by  one  of  the  ordered  lifetimes 

and  approximating  the  conditional  mean  life,  uCt...) 
t  . 

=  r  ^^^R(x)dx,  by  the  total  time  on  test  statistic, 

0*^ 

T(t(i)),  divided  by  n.  Theoretically,  for  infinite  tailed 
distributions,  a  plot  of  the  conditional  mean  life  should 
increase  asymptotically  to  the  mean  life.  Thus,  (from 
the  geometry  in  Figure  10)  as  long  as  the  replacement 
time,  K,  is  finite,  the  optimal  Blind  Replacement  interval 
will  be  finite,  i.e.  in  Blind  Replacement,  we  always 
replace  eventually  regardless  of  the  underlying  distribu¬ 
tion.  So,  if  this  graphical  technique  results  in  an  esti¬ 
mate  tjjj  “^(n)  then  we  replace  at  t^^^j  rather  than 

forego  preventive  replacement  as  we  would  in  Age  Replace- 

/ 

roent. 

To  evaluate  possible  bias  in  our  estimator,  we 
simulated  1000  iterations  of  this  graphical  solution  tech¬ 
nique  for  various  sample  sizes  from  selected  IFR  and  DFR 
Weibull  distributions  and  a  uniform  distribution.  Our 
results  indicate  that  for  a  sufficiently  small  replacement 
time,  K,  =t*  will  converge  to  t*  from  above  as  the 

sample  size  increases,  i.e.,  t*<  t*.  For  a  very  large  K, 
however,  convergence  may  be  from  below.  The  percent  rela¬ 
tive  error  in  the  objective  function,  G(t*),  due  to  the 
biased  estimate  is  provided  below.  For  the  range  of  cases 
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considered,  this  percent  relative  error  is  within  10 
percent. 

For  an  exponential  distribution,  we  can  note  in 
Table  11,  convergence  from  above  for  replacement  times  as 
large  as  three  times  the  mean  lifetime.  The  rapidity  of 
convergence  appears  to  be  a  function  of  the  replacement 
time,  with  convergence  to  within  5  percent  of  the  actual 
optimal  interval  for  sample  sizes  as  small  as  30.  As 
indicated  in  Tables  12  and  13,  convergence  for  IFR  dis¬ 
tributions  may  be  more  rapid  than  for  DFR  distributions 
with  the  Scune  mean  lifetime.  For  these  distributions 
there  is  a  concave  ordering  such  that  conditional  mean 
life  plot  of  IFR  distributions  are  more  concave  than  con¬ 
ditional  mean  life  plots  of  DPR  distributions.  Thus,  the 
geometry  of  the  direct  numeration  method  leads  to  rapid 
converge/ice  for  increasingly  IFR  distributions.  Results 
for  the  uniform  distribution,  shown  in  Table  14,  also 
demonstrate  convergence  from  above  for  small  K  and  conver¬ 
gence  from  below  for  large  K. 

Of  interest  is  the  possibility  of  mitigating  the 
bias  in  this  direct  numeration  technique.  For  a  small 
sample  (say  n  less  than  20)  and  a  small-to-moderate 
replacement  time,  K,  we  might  presume  that  our  estimate 
might  be  greater  them  t*.  Of  course,  if  we  collect  more 

A 

data,  our  estimate,  t*,  of  the  optimal  interval  would  con¬ 
verge,  i.e.  there  would  be  less  bias  in  the  estimate. 
However,  without  any  additional  data  collection,  we  could 
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reaccomplish  our  graphical  technique  using  the  survivor 
time  on  test  statistic  divided  by  n  as  an  approximation 
of  the  conditional  mean  life  in  lieu  of  the  total  time  on 
test  statistic. 

Direct  Numeration  with  Survivor 
Time  on  Test  Statistic 

We  may  be  able  to  improve  on  the  estimate  of  t* 
by  using  the  survivor  time  on  test  statistic,  discussed  in 
Chapter  I,  in  conjunction  with  the  total  time  on  test  sta¬ 
tistic.  Our  strategy  is  to  approximate  the  conditional 

mean  life  by  the  survivor  time  on  test  statistic  divided 

^  i 

by  n  and  identify  as  t  the  point  of  tangency  correspond¬ 
ing  to  the  tangent  line  from  the  point  (-K,0) .  One  naive 
approach  is  to  estimate  t*  as  the  average  between  t  and 

A  A  I 

t*.  For  some  small  replacement  times,  t  may  not  converge 
to  the  optimal  replacement  interval  from  below  because  of 
the  geometry  of  the  direct  numeration  method  (see 
Table  15) .  For  large  n,  of  course,  there  is  no  eppreci- 
able  difference  between  t"  and  t*. 

The  composite  estimator  t  =  (t* +  t  )/2  can  produce 
a  smaller  relative  error  than  the  estimator  t*  as  we  show 
in  Table  16  for  the  exponential  distribution.  This  advan- 

A 

tage  of  the  composite  estimator,  t^,  increases  as  the 
replacement  time,  K,  increases;  however,  it  diminishes  as 
the  sample  size  increases  because  of  the  convergence  of 
t*  and  t  .  A  prudent  strategy,  then,  would  be  to  use  such 
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TABLE  15 

EXPECTED  VALUE  OF  THE  OPTIMAL  BLIND  REPLACEMENT  ESTIMATOR, 
E{T*},  FOR  OBSERVATIONS  FROM  AN  EXPONENTIAL  DISTRIBUTION 
(MEAN  =  200)  WITH  DIRECT  NUMERATION  USING  THE 
SURVIVOR  TIME  ON  TEST  STATISTIC 


E{T*}  for 

Replacement  Time,  K 

= 

Size 

10 

50 

100 

200 

400 

10 

63.4 

122.6 

160.5 

208.4 

261.4 

20 

61.8 

124.4 

167.0 

219.4 

281.6 

30 

61.4 

125.2 

168.0 

222.3 

287.8 

40 

60.9 

126.2 

169.8 

225.1 

291.2 

50 

60.9 

126.7 

170.1 

226.0 

292.9 

60  / 

60.7 

126.0 

170.5 

226.6 

294.7 

70 

60.5 

125.9 

169.9 

226.2 

294.8 

80 

60.5 

125.9 

170.5 

226.4 

295.5 

90 

60.5 

126.1 

170.3 

226.8 

296.0 

100 

60.2 

126.2 

170.7 

227.1 

297.1 

CD 

60.1 

126.5 

171.6 

229.3 

301.1 

TABLE  16 


rM  >  p- 
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a  composite  estimator  regardless  of  the  sample  size  or  the 
replacement  time. 

First  Order  Solutions 


Another  geometric  approach  to  solving  Blind 
Replacement  uses  the  concave  properties  of  G(t)  and  associ¬ 
ated  first-order  sufficient  conditions  expressed  in  equa¬ 


tion  (4.2).  We  can  overlay,  on  a  plot  of  the  conditional 
t 

mean  life,  y=  f  R(x)dx,  a  plot  of  the  function 
(T 

y=  (t+K){R(t)}.  The  latter  function  is  the  left  hand  side 
of  equation  (4.2)  while  the  former  is  the  right  hand  side 
of  the  equality.  The  solution  to  equation  (4.2),  then, 
is  that  t*  where  these  two  plots  cross  as  in  Figure  11. 
Radner  and  Jorgenson  (1962)  show  that  equation  (4.2)  has  a 
vinique  solution.  Therefore,  there  can  be  only  one  cross¬ 
ing  of  these  two  plots.  Since  at  t  =  0,  the  left  hand  side 
of  equation  (4.2)  is  greater  than  the  right  hand  side  (as 
long  as  K  >  0) ,  it  follows  that  the  conditional  mean  life 
plot  will  cross  the  plot  y=  (t+K){R(t)}  from  below.  Given 


.N 


a  set  of  ordered  lifetime  data,  ,  ...,  we  can 

approximate  both  of  these  plots  in  Figure  11  with  lines 
connecting  each  function  evaluated  at  ,  i=l,  ...»  n. 


I 


First,  we  approximate  the  conditional  mean  life, 
t  . . 

f  '^'R(x)dx,  by  the  total  time  on  test  statistic  divided 
0*' 

by  n.  Then,  the  second  function  can  be  approximated  as 
y=  (t  +K)  I  (n-i+1) /n)  using  the  California  plotting  posi- 
tion.  Thus,  we  can  estimate  t*  by  t®  where  t'  is  that 


Pig.  11.  Graphical  First-Order  Solution 
to  Blind  Replacement 


t^j  closest  to  the  crossing  of  the  piece-wise  linear 
approximations  of  both  plots. 

Since  we  approximate  two  functions  with  the  first- 
order  solution  technique  rather  than  one  function  as  with 
a  direct  numeration  solution,  the  first-order  methodology 
is  more  unstable  than  direct  numeration.  This  instability 
is  apparent  in  a  comparison  of  a  first-order  solution 
versus  direct  numeration  as  shown  in  Table  17.  Accord¬ 
ingly,  we  do  not  recommend  using  a  first-order  solution 
technique. 


Summary 

The  Blind  Replacement  model  espoused  by  Radner  and 
Jorgenson  is  a  simple,  probabilistic  representation  of  the 
maintenance  of  a  single-part  item  where  the  objective  is  to 
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TABLE  17 

COMPARISON  OF  EXPECTED  VALUES  OF  OPTIMAL  BLIND  REPLACEMENT 
ESTIMATORS  FOR  AN  EXPONENTIAL  DISTRIBUTION  (MEAN  =  200) 
USING  A  DIRECT  NUMERATION  (t*)  TECHNIQUE  VERSUS  A 
FIRST  ORDER  (T§)  SOLUTION  TECHNIQUE 
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maximize  the  item*  s  availability  when  the  item' s  condi¬ 
tion  is  unknown.  This  model  is  eunenable  to  graphical 
solutions  similar  in  nature  to  Bergman's  technique  for 
Age  Replacement.  He  proposed,  herein,  several  solution 
techniques  and  discussed  the  rapidity  of  convergence  of 
each  to  the  optimal  Blind  Replacement  interval  based  on 
simulations  from  several  common  life  distributions.  He 
note  that  the  proposed  first-order  solution  is  relatively 
unstable  compared  to  direct  numeration  and,  therefore, 
do  not  recommend  its  use.  As  for  direct  numeration  solu¬ 
tion  techniques,  we  show  that  the  use  of  the  total  time  on 
test  statistic  leads  to  biased  estimates  of  t*.  This  bias 
might  be  reduced,  for  moderate  replacement  times,  by  using 
a  combination  of  graphical  estimators  based  on  the  total 
time  on  test  statistic  and  the  survivor  time  on  test  sta¬ 
tistic. 


CHAPTER  V 

CONFIDENCE  BOUNDS  IN  MAINTENANCE  PLANNING 

Given  a  Scunple  of  failure  data,  the  graphical  tech¬ 
niques  for  Age,  Block,  and  Blind  Replacement,  discussed  in 
Chapters  II,  III,  and  IV,  provide  point  estimates  of  the 
optimal  replacement  interval-  Ideally,  we  would  like  to 
see  confidence  bounds  about  these  point  estimates  in  order 
to  judge  the  degree  of  dispersion  inherent  in  any  sampling 
technique.  However,  for  a  single  sample  of  size  n,  the 
distribution  of  the  optimal  replacement  interval  estimate, 
t^jj  ~  t*,  is  that  of  the  order  statistic  where  J  is 
itself  a  random  variable.  Its  sampling  distribution  is  of 
an  unknown  form  and  is  not  amenable  to  straightforward 
statements  of  statistical  variability.  Thus,  measures  of 
dispersion,  which  gage  sampling  accuracy  for  small  scunple 
sizes,  may  have  to  be  obtained  by  artificial  sampling 
experiments.  In  this  chapter,  we  discuss  several  tech¬ 
niques  for  bounding  the  true  optimal  replacement  interval, 
t* ,  using  a  distribution  free  resampling  concept  known  as 
the  "percentile  method."  We  illustrate  this  methodology 
for  Age  Replacement  and  Blind  Replacement;  however,  this 
resampling  concept  is  also  appropriate  for  Block  Replace¬ 
ment. 


Efron  (1982)  discusses  a  resampling  procedure 
known  as  the  "bootstrap"  and  advocates  the  percentile 
method,  a  nonpar aunetric  confidence  interval  technique 


developed  from  the  bootstrap  procedure.  The  bootstrap  is 
a  computational  process  for  estimating  attributes  of  some 
general  random  variable,  X,  with  an  unspecified  cumulative 
distribution  function,  F(x).  This  process  involves 
repeated  resampling  of  a  single  set  of  observations, 

{Xj^,  x^,  ...,  where  n  resample  observations  are 

uniformly  drawn  with  replacement  from  the  set  of  the 
original  n  observations.  Each  bootstrap  resample  set, 

{x^,  x^,  ...,  can  be  used  to  compute  some  aspect  of 

X,  while  repeated  bootstrap  resampling  provides  statistics 
on  this  aspect.  In  our  case,  we  want  to  estimate  t*  from 
a  set  of  n  observed  lifetimes,  {t,,v<»»*<t,  and 
specify  a  confidence  bound  about  t*  without  collecting 
more  data.  Accordingly,  we  can  resample  with  replacement 
from  the  original  lifetime  data  and  accumulate  b  sets  of 
bootstrap  data.  It  j  <  . . .  <  t  }  ...  <  . . .  <  t'^j  }. 

For  each  of  the  b  sets  of  bootstrap  observations,  we  can 
ccxnpute  t*^  (i=l,  ...,  b) ,  the  i*"^  bootstrap  estimator  of 
t*.  The  collection  of  b  bootstrap  estimators  can  them¬ 
selves  be  ordered  as  { t*  £  . . .  £  t*  ,  and  a  100  (1-a) 
percent  confidence  bound  can  be  constructed  frcxii  these 
ordered  bootstrap  estimators  as  t*  ] . 

We  refer  to  this  confidence  bound  as  the  "percentile 
method"  bound. 


In  order  to  investigate  the  performance  of  the 
percentile  method  for  Age  Replacement  and  Blind  Replace¬ 
ment,  we  turn  to  Monte  Carlo  simulation  involving  several 
distributions  from  the  Weibull  family  whose  true  optimal 
replacement  intervals,  t*,  are  known  a  priori. 


The  Percentile  Method  for 


Age  Replacement 

We  applied  Efron's  percentile  method  for  confi¬ 
dence  bounds  to  Bergman's  graphical  Age  Replacement  esti- 
mator,  t*,  based  on  50  observations  from  selected  IFR 
Weibull  distributions,  each  with  the  same  mean  value. 

Our  analysis  is  based  on  1000  iterations  of  50  original 
observations  per  iteration  from  the  respective  Weibull  dis¬ 
tributions.  At  each  iteration,  we  resampled  the  50  obser¬ 
vations  (an  arbitrary)  100  times  with  a  bootstrap  sample 

size  of  50.  If  the  last  ordered  observation  in  the  boot- 

) 

strap  sample  minimized  the  objective  function,  then  for 
numerical  reasons  we  estimated  t*  as  a  large  number  in 
lieu  of  using  infinity.  We  then  constructed  various  con¬ 
fidence  level  bounds  for  each  iteration,  and  noted  whether 
the  bound  contained  the  true  optimal  Age  Replacement  inter¬ 
val,  t*,  an  appropriate  fraction  of  the  time.  The  results 
from  this  simulation,  summarized  in  Tables  18  through  20, 
indicate  that,  for  some  distributions  and  cost  ratios,  the 
percentile  method  may  not  provide  confidence  bounds  that 
are  wide  enough  to  contain  t*  with  the  correct  probability. 
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When  there  are  large  differences  between  the  cost 
of  failure  and  the  cost  of  replacement  (i.e  a  small  cost 
ratio) ,  the  percentile  method  produces  confidence  bounds 
which  have  a  low  probability  of  containing  t*  and  the 
potential  for  large  relative  error  in  the  objective  func¬ 
tion.  A  partial  explanation  of  this  result  involves  the 
geometry  of  Bergman's  technique  which  leads  to  the  first- 
order  statistic  as  an  estimate  of  t*  when  the  cost  ratio 
is  small.  With  scant  original  samples,  it  is  possible 
that  the  smallest  observed  failure  time  is  greater  than  t* 
thereby  leading  to  an  inability  of  the  percentile  method 
to  capture  t*.  This  difficulty  can  be  diminished  by 
increasing  the  original  sample  size.  Thus,  for  small  cost 
ratios,  a  prudent  course  of  action  would  be  to  collect 
reasonably  large  numbers  of  observed  failure  times. 

yfith  moderate  cost  ratios  (i.e.  a  moderate  differ¬ 
ence  between  the  cost  of  failure  and  the  cost  of  replace¬ 
ment)  ,  the  percentile  method  performs  adequately  for 
Weibull  distributions  with  shape  parameters  greater  than 
or  equal  to  2.0.  For  large  cost  ratios  (which  implies  a 
nominal  difference  between  the  cost  of  failure  and  the  cost 
of  replacement)  ,  the  percentile  method  wor]cs  well  as  long 
as  the  Weibull  distribution  is  not  too  IFR.  Although  the 
90  percent  confidence  intervals  for  these  moderate  to  large 
cost  ratios  can  be  very  wide  (as  shown  in  Table  19) ,  the 
average  maximum  relative  error  in  the  objective  function 
within  these  intervals  is  less  than  10  percent. 


The  application  of  the  percentile  method  to  Age 
Replacement  can  result  in  a  confidence  bound  that  is  nil 
The  likelihood  of  this  possibility  increases  as  the  sta¬ 


tistical  confidence  of  the  bound  decreases.  For  the 
reasonably  large  confidence  bounds  (greater  than  70  per¬ 
cent)  ,  we  never  encountered  a  zero  width  confidence  bound 
throughout  our  simulation  activity.  Thus,  for  sufficiently 
large  numbers  of  observations,  say  at  least  50,  there 
appears  to  be  minimal  risk  of  a  zero  range  confidence 
bound  as  long  as  we  consider  confidence  levels  greater  or 
equal  to  70  percent. 

The  Percentile  Method  for 
Blind  Replacement 

We  also  applied  Efron's  percentile  method  for  con¬ 
fidence  bounds  to  our  TTT  statistic  direct  numeration 

A 

estimator,  t*,  based  on  50  observations  from  selected  ZFR 

f 

Welbull  distributions,  each  with  the  same  mean  value. 

Again  our  analysis  is  based  on  1000  iterations  of  50 
original  observations  from  the  respective  Welbull  dis¬ 
tributions.  At  each  iteration,  we  resampled  the  50  obser¬ 
vations  100  times  with  a  bootstrap  sample  size  of  50,  con¬ 
structed  various  confidence  bounds  for  each  Iteration,  and 
noted  whether  the  bound  contained  the  true  optimal  Blind 
Replacement  interval.  The  results  from  this  simulation, 
contained  in  Table  21,  indicate  that  the  percentile  method 
produces  confidence  bounds  which  are  not  quite  wide  enough 
to  contain  t*  with  the  correct  probability. 
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The  reason  for  this  weak  performance  involves  the 
positive  bias  in  the  Blind  Replacement  estimator  t*  (using 
the  total  time  on  test  statistic)  which  we  noted  in 
Chapter  IV.  One  approach  for  improving  this  performance 
is  to  modify  the  percentile  method  by  reducing  the  boot¬ 
strap  resample  size  to  some  fraction  of  the  original 
observed  sample  size.  Decreasing  the  bootstrap  resample 
size  will  increase  the  variability  in  our  biased  estimator 
and,  ideally,  will  widen  the  confidence  bound.  Unfor¬ 
tunately,  this  approach  does  not  work  for  Age  Replacement 
as  it  does  for  Blind  Replacement. 

We  replicated  our  Blind  Replacement  simulations 
with  a  modified  percentile  method  involving  bootstrap 
resample  sizes  of  80  percent  of  the  original  sample  size 
in  lieu  of  the  100  percent  bootstrap  resample  size  used 
above.  In  other  words,  each  of  the  100  bootstrap  resampled 
data  sets  contained  only  40  observations  drawn  with 
replacement  from  the  original  50  observations.  The 
results,  depicted  in  Tables  22  through  24,  are  quite 
satisfactory.  This  modified  percentile  method  produced 
confidence  bounds  which  contained  t*  approximately  the 
correct  number  of  times  in  1000  iterations  as  can  be  seen 
in  Table  22.  The  90  percent  confidence  bounds  given  in 
Table  23  are  sufficiently  narrow  to  be  able  to  bracket 
t*  to  within  a  10  percent  maximum  relative  error  in  the 
objective  function  as  shown  in  Table  24. 
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Nonparametric  confidence  bounds  for  the  optioial 
replacement  interval,  t*,  can  be  constructed  via  a 
resan^ling  technique  which  Efron  (1982)  calls  the  per¬ 
centile  method.  This  method  uses  the  order  statistics 
generated  from  a  bootstrap  resaiiq>ling  scheme  to  form 
bounds  with  a  given  level  of  statistical  confidence.  A 
Monte  Carlo  simulation  of  graphical  solution  techniques 
for  Age  Replacement  and  Blind  Replacement  suggests  that 
in  many  cases  the  percentile  method  may  not  produce  bounds 
with  the  correct  statistical  confidence.  Moreover,  in 
some  cases,  the  average  maximum  error  in  the  objective 
function  within  the  confidence  bounds  can  be  substantial. 
For  Blind  Replacement,  ho%#ever,  a  modified  percentile 
method  can  produce  bounds  with  satisfactory  levels  of  sta- 
tistical/conf idence  and  nominal  relative  errors  in  the 
objective  function. 


APPENDIX 


LIST  OF  SPECIAL  NOTATION 


E{V} 

Expected  cost  during  a  typical  renewal 
cycle  (p.  4) 

E{L} 

Expected  renewal  cycle  length  (p.  4) 

F(t) 

Failure  distribution  (p.  5) 

R(t) 

Survival  distribution  (p.  5) 

M(t) 

Conditional  mean  life  (conditioned  on 
censoring  lifetimes  at  time  t) (p.  5) 

t* 

Optimal  replacement  interval  (p.  6) 

Ordered  observations  where 

T(t^i)) 

Total  time  on  test  statistic  (p.  7) 

Scaled  total  time  on  test  statistic 
(p.  8) 

S(t(i,) 

Survivor  time  on  test  statistic  (p.  15) 

Sample  path  of  a  renewal  process  (p.  16) 

M(t) 

Renewal  function  (p.  16) 

t* 

Estimate  of  t*  (p.  31) 

E{T(i)} 

Expected  value  of  the  i^^  order  sta¬ 
tistic  (p.  32) 

E{T*} 

Expected  value  of  the  estimator  of  t* 

(p.  39) 

m(t) 

Renewal  density  function  (p.  43) 

Estimate  of  t*  based  on  i^^  bootstrap 
resample  of  n  observed  lifetimes  (p.  75) 
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